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■ The formation of superheavy nuclei in "^^Ca+^^^T/i, ^^*[/, 242,244p^ '^^^Cm 

reactions and their subsequent decay are studied within the quantum mechani- 
cal fragmentation theory (QMFT) and the QMFT based preformed cluster-decay 
model (PCM) of Gupta and collaborators. According to QMFT, all these ^*Ca- 
induced reactions are cold fusion reactions with relative excitation energies larger 
than for the P6- induced cold fusion reactions and smaller than for the lighter beam 
^.^ , i.e. Mg, Si or S-induced hot fusion reactions. The same reactions were first sug- 

! gested by Gupta et al. in 1977 on the basis of QMFT, and this study re-establishes 

Q\ ■ the same result. In fact, for such heavy isotopes of Z=110 to 116, ^^Ca is shown to 

be a better beam for cold fusion, but ^°Ca is a radioactive nucleus. The a-decay 
half-lives of these nuclei after 3n and/ or 4n evaporations, i.e. of the evaporation 
residues of these compound systems, calculated on PCM compare reasonably well 
with experiments published by Dubna group and another recent calculation. As 
^ ' expected for such rare decays, PCM calculations show that the a-preformation 

I factors are small, ~ 10^^ to 10^^*^. The possible competition of a-decays with 

heavy cluster emissions from these superheavy nuclei is also probed from the point 
of view of searching for new nuclear structure information and possible future ex- 
periments with such exotic nuclei. The decay half-lives for some clusters are in fact 
shown to be lower than the limits of experiments for nuclei with enough available 
atoms. 
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1 Introduction 



The formation of new and superheavy elements through the use of highly neutron 
rich beam of ^^Ca on neutron rich actinide targets, such as '^'^^Pu, '^'^^Cm and 
'^^'^Cf, was first suggested by Flerov [1] and was used at his JINR Dubna U- 
300 heavy-ion cyclotron for the successful synthesis of '^^'^No in a cold (17-18 MeV 
excitation) 2n emission reaction with ^°^P6 target [2] . However, the same could not 
be continued to heavier elements due to the poor sensitivity level of the equipment 
at that time [3] . Theoretically, Nix [4] later suggested the use of very asymmetric 
target-projectile combinations, such as ^^°Cm and ^^'^Fm with ^^Ca, on the basis 
of the kinetic energy calculations made for different nuclear shape configurations by 
using an idealised liquid drop model. In a more complete theory, called quantum 
mechanical fragmentation theory (QMFT) [5, 6, 7, 8, 9, 10, 11, 12], based on two 
centre shell model and Strutinsky renormalization procedure, Gupta, Sandulescu 
and Greiner [11] also suggested the use of ^^Ca beam, on the minimum energy 
considerations, with even-even 19^-204^^^ 204-208p^^ 232^^^ 234-238^;^ 24o-244p^ 

and ^^^"^"^^Cm targets for synthesizing the various isotopes of even Z=100, 102, 
and 110-116 (only even Z nuclei were studied). The calculations also showed that 
for the heavier mass isotopes at the minimum in energy either the projectile is 
not ^^Ca or the corresponding target nucleus is not stable in nature. On the 
other hand for more neutron-deficient isotopes, the potential energy surfaces tend 
to become flat and ^^Ca is no more a favoured nucleus since other new minima 
start to develop. Furthermore, for Z=104-108, the targets to be used with ^^Ca 
projectile are the unstable, very difficult to handle Po, Rn and Ra nuclei. The 
key consideration behind the QMFT is the shell closure effects of one or both 
the reaction partners and hence of "cold reaction valleys" or the reaction valleys 
leading to cold fusion. Four or five such cold reaction valleys were always found 
to exist, '^^Ca being one of them. Their relative excitations would mean the cases 
of "cold", "warm/tepid" and "hot" fusion reactions. In other words, all cases 
of potential energy minima are in fact of "cold" reactions with different relative 
excitation energies, and the real "hot" reactions are the ones coming from outside 
the potential energy minima. 

More recently, an intense ^^Ca beam at low consumption of material in the 
ion source was developed [13] and used with the upgraded Dubna cyclotron U- 
400 and new recoil mass separators (VASSILISSA or GNS) for the formation of 
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superheavy elements 110 to 116 [14, 15, 16, 17, 18, 19]. The targets used are 
^^'^Th, "^^^U, 242,244p^ 24:8(j^^ ^^^^^ f^j. ^^^j. ^]^g Coulomb barrier energies, the 

excitation energy E* ~30-35 MeV, in between the one for cold {E* ~ 10-20 MeV) 
and hot {E* ~40-50 MeV) fusion reactions. Thus, the ^^Ca-induced reactions are 
referred to as "warm/ tepid" fusion reactions, in between the "cold" and "hot" 
fusion reactions. For some ^^Ca-induced reactions, E* ~40 MeV also. Thus, the 
resulting compound systems de-exite by 3n and/ or 4n evaporations, compared to 
In and 2n in cold and 5n in hot fusion reactions. The cold fusion reactions are 
based on Pb and Bi targets with Ti to Zn beams, and the hot fusion reactions 
use heavy actinide targets from Th to Cf with very light C to 5" projectiles. In 
all the "warm" fusion reactions studied here, the final nuclei formed after the 
evaporation of neutrons (and 7-ray emission), the evaporation residues EVRs, are 
relatively long-lived nuclei and decay only via a-particles, ending in spontaneous 
fission (SF) of the last nucleus at the edge of the stabihty region. This is the 
complete fusion reaction channel, giving the a- genetically related nuclei, called 
a-decay chain. Another decay channel for the excited compound system is that 
of fission-like decays: the fusion-fission, quasi- fission, etc., which has also been of 
interest both experimentally [20, 21, 22, 23] and theoretically [24, 25, 26]. It may 
be relevant to mention here that the fusion-fission process in hot superheavy nuclei 
could also be treated as a dynamical cluster decay process [27] advanced recently 
for light and medium mass hot compound systems. 

In this paper, we first re-investigate the use of the ^^Ca beam for the synthesis 
of superheavy elements within the dynamical fragmentation theory via cold fusion 
reactions. This has become essential because the binding energy data has improved 
considerably since its last use [11] by some of us in 1977. Specifically, we consider 
the compound systems ^^"llO*, '^^^112*, '^^0,292^^* '^^f^UQ*^ formed in reactions 
of ^«Ca beam on ^^^Th, ^^^U, 242,244p^ 248^^ targets at various incident 
energies. As already stated above, the excitation energy E* — 32-35 MeV in 
these reactions and hence they all de-excite by 3n and/ or 4n evaporations, to 
give 277110, 283112, 287,288,289114 292^^^ parents. Then, as a second aim of this 
paper, we investigate the observed a-decay characteristics of these nuclei within the 
preformed cluster decay model (PCM) of Gupta and Collaborators [28, 29, 30, 31, 
32, 33] which is also based on the dynamical fragmentation theory. Furthermore, a 
possible branching of a-decay to some (theoretically) most probable heavy cluster 
decays is also studied with a view to look for some new or known magic daughters 
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in this exotic process of, so-called, cluster radioactivity. 

The paper is organised as follows. A brief description of the fragmentation 
theory and PCM is given in section 2. The results of our calculations, compared 
with Dubna experiments and other recent works, are presented in section 3. A 
summary of our results is added in section 4. 

2 The Theory 

2.1 Quantum mechanical fragmentation theory 

The quantum mechanical fragmentation theory (QMFT) is a dynamical theory 
of heavy-ion collisions, specifically the three cold processes of fusion, fission and 
cluster radioactivity, worked out in terms of the coordinates of mass (and charge) 
asymmetry r] = (Ai — A2)/A (and rjz = {Zi — Z2)/Z), the relative separation 
distance R, the deformations /3i and of two nuclei (or, in general, fragments), 
and the neck parameter e. Taking motions in f] and r]z as weakly coupled, the 
time-dependent Schrodinger equation in rj, 

HmM = in^^>M, (1) 

is solved for R(t) treated classically and the other coordinates ^i, ^2 and e fixed by 
minimizing the collective potential in these coordinates [34, 35] . We find that for a 
given compound system (formed by different target -|- projectile combinations), a 
few nucleon to a large mass transfer occurs for the target -|- projectile combinations 
coming from outside the potential energy minima (the real "hot" combinations), 
whereas the same is zero (no transfer at all) for the target -|- projectile combinations 
referring to potential energy minima (the "cold" combinations); see Fig. 14 in [35]. 
This means that for cold reaction partners, the two nuclei stick together and form a 
deformed compound system, as is illustrated in Fig. 2 of [9]. On the other hand, a 
few nucleon transfer may occur if only a "conditional" saddle is formed (see Figs. 1 
and 2 in [36]). Eq. (1) is solved for only a small number of heavy systems, since its 
solution is very much computer-time consuming. However, certain simplifications, 
discussed in the following, have resulted from actual calculations performed over 
the years which seem to work rather nicely. 

The dynamical QMFT estabhshes the best target -|- projectile combination, 
with the variation of relative separation coordinate R allowing to pass continuously 
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from a separated pair of nuclei to a cool compound nucleus [8, 9] or any other 
process like fission or deep inelastic collision, etc., [25, 26]. However, the potentials 
V(R,r^) and V(R,?72), calculated within the Strutinsky renormalization procedure 
iy—^LDM + (^U) by using the appropriate liquid drop model for Vldm [37] and 
the asymmetric two-center shell model (ATCSM) of Maruhn and Greiner [38] for 
the shell effects 5\J, show that the motions in both r] and rjz are much faster than 
the R-motion. This means that both the potentials V(R,r/) and V(R,r;z) sue nearly 
independent of the R-coordinate (see e.g. Fig. 1 in [8] or [24]) and hence R can be 
taken as a time-independent parameter. Then, the time-dependent Schrodinger 
equation (1) in rj reduces to the stationary Schrodinger equation in rj, 

where R is fixed at the post-saddle point, justified by many good fits to both 
fission and heavy-ion collision data [39] and by an explicit, analytical solution of 
time-dependent Schrodinger equation in r)z coordinate [40]. An interesting result 
of these calculations is that the mass (and charge) distribution yields, 

y(A) =1 Mv{A^)) r (i=l,2) (3) 

(and Y{r)z)) are nearly insensitive to the detailed structure of the cranking masses 
B^^, calculated consistently by using ATCSM. The nuclear temperature effects are 
included here through a Boltzmann-like function 

oo 

I r= E I r expi-E^^^/n (4) 

i/=0 

with the temperature T (in MeV) defined, via the compound nucleus excitation 
energy, as 

E*^E,m + Qin^^T'-T. (5) 

Qin is the Q-value for reaction partners. Also, the shell corrections SU are T- 
dependent, but this could not be included here since we are using the experimental 
binding energies where SU are there in them but could not be separated out in a 
model independent way. However, the T-values involved here are small, ~1 MeV 
only. 

Thus, the static potentials V(?7) (and V(?7z)) contain all the important infor- 
mation of a colliding or fissioning system. Furthermore, since these potentials are 
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nearly independent of the choice of R- value (as discussed above) , they are calcu- 
lated at some critical distance R — Ct — Ci -\- C2 where the two nuclei come in 
close contact with each other. This means fixing the neck parameter e—1 and the 
potential V(?7, r]z) given simply as the sum of the ground state binding energies 
Bi of two nuclei, the Coulomb interaction E^. between them plus the additional 
attraction due to nuclear proximity Vp [41], 

y (a, 77, vz)^-jz B{A, z,, A) + + Vp. (6) 
i=i 

The Cj, defining Ct, are the Siissmann central radii Cj = Ri — {1/Ri), with the radii 
Ri = 1.28^-/^ - 0.76 + 0.8^"^/ V"^- The binding energies Bi are from the 1995 
experimental compilation of Audi and Wapstra [42] and from 1995 calculations 
of MoUer et al. [43] whenever not available in [42]. The charges Zi and Z2 are 
determined by minimizing the potential in r]z coordinate, which automatically 
minimizes the ^ coordinates. Note that the minimized ^'s are not always for the 
spherical nuclei since the total binding energy Bi + B2 of the reaction partners 
is minimized and not their individual Bi or B2. Also, it may be pointed out 
that Vp, the additional attraction between nuclear surfaces, was not added in our 
earlier calculations, whose effect is to change the relative heights of the potential 
energy minima, and hence the relative excitations, but not their positions [36] . The 
positions of the minima are due to shell effects only. Also, the role of deformation 
in both Ec and Vp is neglected here because their combined effect is shown [44] 
to lower the interaction barriers but not the relative formation yields. In other 
words, Eq. (6), without Vp, formed the basis of our first calculation on "cold 
fusion" reaction valleys [8, 10, 11, 12], which was later optimized by adding the 
requirements of smallest interaction barrier, largest interaction radius and non- 
necked (no saddle) nuclear shapes [9]. Like necked- in shapes are known [6, 7] to 
witness the preformation of fission fragments, non-necked shapes are the signatures 
of cold fusion of two nuclei [9] . The contribution due to Vp is, however, added in 
the present calculations. 

2.2 Preformed cluster-decay model 

The preformed cluster-decay model (PCM) is based on the QMFT and hence uses 
the same coordinates as are introduced above. In a PCM, the decay constant A, 
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or the decay half-life T1/2, is defined as 

A = 1^ = Po^oP, (7) 

-'1/2 

with Pq as the cluster (and daughter) preformation probability in the ground state 
of nucleus and P the barrier penetrability which refer, respectively, to the 77 and 
R motions. The uq is the barrier assault frequency. In principle, the coordinates 
77 and R are coupled, but Eq. (7) is always written in a decoupled approximation, 
which further justifies our simplifying conditions of the last subsection. 

The Pq are the solutions (3) of the stationary Schrodinger equation (2) in 77 
for the ground-state u — 0, i.e. Pq — Y°{r]). The mass parameters B^^{ri) are the 
classical hydrodynamical masses of Kroger and Scheid [45] since the two masses 
(cranking and hydrodynamical) are shown to give similar results for heavier nuclei. 

The P is the WKB tunneUing probabihty, calculated for the tunnelling path 
shown in Fig. 1, as P = PjPb with 

Pi = exp[-l r {2i,[V{R) - V{Ri)]y/'dR] (8) 

Pb = exp[-l r'{2/i[y(P) - Q]y/'dR], (9) 

with the first turning point Ra = Ct and the second turning point Rb defined by 
V{Ri,) — Q-value for the ground-state {a or cluster) decay. These integrals are 
solved analytically [29]. This choice of Ra — Ct{= Ci + C2), instead of P = Pq, 
the compound nucleus radius, assimilates to a good extent the effects of both 
the deformations of two fragments and neck formation between them [44]. In 
other words, the two-centre nuclear shape is simulated here through a neck-length 
parameter which for actinides is nearly zero [44] . We have taken it to be the same 
for superheavy nuclei. The role of deformation in V(R) is to lower the interaction 
barriers [44], which is achieved here by raising P^ from Pq to Cf. 
The assault frequency uq in (7) is given simply as 

(10) 

with E2 — {Ai/A)Q, the kinetic energy of the hghter fragment, for the Q- value 
shared between the two products as inverse of their masses, /i is the reduced mass 
here. Eq. (10) usually results in uq ?=:i 2.7 x lO^^s"^, whereas the more often used 
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value in literature is higher for even parents and lower for odd parents [32]. We 
choose to use here uq = 2.7 x lO^^s"^ for odd parents and uq — 2.7 x lO^^s"^ for 
even parents. 

3 Calculations 

The superheavy nucleus formation process: Figure 2 shows our calculated fragmen- 
tation potentials V{rj) at a fixed R — Cj-value, based on Eq. (6), for some excited 
compound systems ^^^16*, "^^HU*, "^^^12* and 280^10* formed via ^^Ca-induced 
reactions. We find in Fig. 2 that in each case there is a minimum corresponding 
to ^°Ca nucleus (the ^^Ca minimum lies in its immediate neighbourhood but at 
somewhat higher energy), which is in addition to other minima at symmetric (or 
nearly symmetric), the ^^^Pb and ^^Kr or their neighbouring nuclei, respectively, 
with Z = 80 ± 2 and 36 ± 2 (the Pb and Kr minima merge into one minimum 
for Z > 116 elements) and the ones corresponding to the super-asymmetric com- 
binations using heavier actinides with much lighter beams. The ^^Ca minimum is 
deepest only for the hghter isotopes of these nuclei, e.g. for ^^°116, ^^^114, ^'''^112 
and ^™110 (see Fig. 1 and Table 1 in [46]). Apparently, all these minima cor- 
respond to at least one closed shell nucleus and hence to "cold" reaction valleys. 
Also, the relative excitation energies are different, the ^^Ca minimum always ly- 
ing higher than the ^°^P6 (and/or ^^Kr) and the symmetric (or near symmetric) 
minima, but deeper than the ones corresponding to many heavier actinides with 
hghter beams of Mg, Si, S, etc.,. The excitation energy can be determined by the 
hight of the minimum with respect to the ground state of the compound system. 
Knowing that temperature effects are not added here in Fig. 2, the approximate 
excitation energy E* ~30 MeV for Ca minima, to be compared with E* ~10 MeV 
for Pb minima and ~40 MeV for heavy actinides. For lighter elements (Z=102 
and 104), the ^^Ca minimum is deeper than the Pb and Kr minima, and its rel- 
ative depth is only 15 to 20 MeV [11, 46]. This means to say that Ca-induced 
reactions for heavy elements {Z > 106) give rise to "cold" compound systems with 
an "intermediate" amount of excitation energy, and hence the name "warm or 
tepid" with respect to "cold and hot" fusion. The number of neutrons emitted in 
^^Ca-induced reactions (3n or 4n) also he in between the cold (In or 2n) and hot 
(5n) fusion reactions. In other words, according to the QMFT, the ^^Ca beams 
are as good as the ^°^P6 beams (or the hghter ones like Mg, Si and S) for forming 
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cold compound nuclei with relatively higher (or lower) excitation energies. The 
calculated interaction barriers and nuclear shapes also support that Ca-induced 
reactions are as good as the P6-induced ones [9, 26, 46]. Finally, it may be relevent 
to note that the ^°Ca nucleus, preferred for cold compound nucleus formation in 
heavier isotopes of superheavy nuclei studied here, is radioactive whereas the ^^Ca 
is not a radioactive nucleus. The predicted use of radioactive nuclei (as a beam 
and/ or a target) for heavier isotopes of superheavy nuclei is a general result of 
the QMFT [46]. 

We have also calculated the fragment mass distribution yields Y{ri), using Eq. 
(3), taking the view that, since fragments related to the minimum in the potential 
Vlrj) are more probable, the yields Y{'r)) must give the intermediate (two) fragment 
formation yields or, in short, the formation yields for a cool compound nucleus [47], 
where the contribution of barrier penetration is not included. Fig. 3 illustrates 
our calculation for two compound systems ^^^116* and ^^^114*. For the mass 
parameters, we have used the classical hydrodynamical masses [45]. Excluding 
the cases of reactions involving symmetric combinations, in view of their forming 
necked-in shapes, the formation yields for ^^Ca-induced (here ^°Ca) cold fusion 
reactions are apparently larger than for the light mass beams but smaller than 
for P6-induced reactions. The contribution of penetrability P would, however, 
quantify these results. 

The a-decay process in superheavy nuclei: Table 1 and Fig. 4 show the results 
of our calculation for a-decay half- lives, in terms of logioT^^^ (s), for the measured 
complete decay chains of ^''"''110, 287,288,289 ^^^^4 "^^^IIQ parents (no a-decay of 
^^^112 parent is observed), compared with the experimental data pubhshed by 
Dubna group and a very recent calculation of Royer and Gherghescu [48, 49] based 
on the generalized hquid drop model (GLDM). The Dubna data is for a-particle 
energies and decay times r, which agree with Qq- values (within less than 
ISOKeV) and a-decay half-lives T^^2 (within less than one order of magnitude), 
respectively, via the well known formula of Viola and Seaborg [50] whose parame- 
ters are fitted to measured Qa and T{^2 values of 58 even-even Z > 82, N > 126 
nuclei [51]. We have, therefore, used this data to compare with the calculated 
Qa and T'{^2"Va'lues. The resulting differences are of 0.2-0.4 atomic mass num- 
bers (compare the resulting 114.4lo;8 ^^id 110.2lo;8 [14] , respectively, with A=114 
and 110). Wherever more than one chain is observed, we present the data for all 
chains. Table 1 also lists the respective Qa- values and the preformation factor Pq 



9 



and penetrability P calculated on PCM. The PCM calculated Q-values {Q'^^'") 
based on experimental [42] and MoUer et al. [43] binding energies, whereas that 
of GLDM on Thomas-Fermi (TF) model [52]. The Q'^^p*- are taken to be the ex- 
perimental £^Q,-values measured in respective a-decay chains. In order to illustrate 
the role of Q- value, we have also calculated the decay half-lives within PCM using 
the experimental Q^^^'^ -values, and compared with the results of our calculations 
using Qa'^''. 

We first notice in Table 1 that a decrease in Q- value by about 0.5 MeV increases 
the Ti-value by about one order. This is also illustrated in Fig. 4 for ^''"''110 
and 288,289-|^24 decay chains. Apparently, in one case (Z=114 decay chains) use of 
experimental Q-values improve the comparisons, whereas in the other case (^^^110 
decay chain) this spoils it. In general, the comparisons of Ti -values for the two 
models with experiments are within experimental errors, i.e. within less than two 
orders of magnitude. In some isolated cases, like for ^^^106 nucleus in ^^^110 decay 
chain, the trends of experiments are different from calculations and hence need 
further checking. In total, the data obtained do not show any strange behaviour 
and hence do not call for any special criticism [53]. 

The interesting result to note in Fig. 4 is that both the GLDM and PCM calcu- 
lations give similar trends, which is due to their used Q-values. This is illustrated 
in Fig. 5, where the calculated Q-values on two models are compared with exper- 
iments. The experimental E'^-values (taken as (5^^^* -values) are in general higher 
and hence are the cause of one to two orders of magnitude deviations between the 
calculated Ti and measured r- values. 

2 

Finally, in Table 1, we have a look at Pq and P, calculated on PCM. The Pq 
are of the order of 10~^ to 10~^° which means that all these a-decays are rare 
in nature, as is actually the case. In GLDM, Pq — 1, which means that the P 
calculated on PCM should be higher, due to either the higher Q-values or lower 
barriers. However, Fig. 5 shows that the Q-values used in PCM are lower than 
those used in GLDM. Also, it may be noted that vq — 10^°s~^ in GLDM, an 
order of one and three lower from that used in PCM, respectively, for odd and 
even parents. All these factors have a combined effect on the nature of predictions 
made by a particular model. 

The heavy cluster- decay process in superheavy nuclei: Figure 6 shows the results 
of our calculations for some heavy cluster decays whose preformation factors Pq 
are the largest and hence refer to the minima in fragmentation potentials V{r]), 
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for each of the parents in the a-decay chain of, say, ^""llO nucleus. The heaviest 
cluster included here in the calculations is with Z=20 (^^~^^Ca) because an earher 
calculation on PCM shows that the two processes of cold fission and cluster decay 
become indistinguishable for clusters heavier than of mass A2 48 [31]. The 
results of other decay chains are not presented hence since this study at present is 
more of an academic interest than for comparisons with experiments. 

First of all we notice in Fig. 6 that, almost independent of the parent mass, 
the Q- value increases as the size of cluster increases. On the other hand, P and Pq 
are further smaller than for a-decay, but present an interesting result. We notice 
that though ^°Be is best preformed in all parents, its P-values are very small. Such 
a result is important since 71/2 is a combined effect of both Pq and P {uq being 
constant). For ^°Be, the Ti/2 is much larger than for other clusters such that the 
studied parents could be said stable against ^°Be decay. In other words, in addition 
to a-decay and fission, ^^C, ^^Si and/ or Ca clusters present the best possible cases 
of cluster decays for any of the parents of ^''"''110 a-decay chain. There seems to 
be no new shell stabilizing effects in the considered decays, except for ^^Be decays 
of ^^^106 or ^^^108 nuclei, possibly due to the known deformed magicity of these 
superheavy nuclei [54]. Interesting enough, the calculated half-lives for ^^~^^Ca 
decays he far below the present limits of experiments, which go upto ~ lO^^s [32] 
for nuclei where enough atoms are available. The closed shell effects of a cluster, 
however, are not yet observed in the exotic cluster radioactivity studies. 

4 Summary of our results 

We have re-investigated the problem of synthesizing superheavy elements via the 
use of ^^Ca beam, within the quantum mechanical fragmentation theory (QMFT) 
that was first used by Gupta, Sandulescu and Greiner [11] in 1977. The most 
probable, specific reaction partners were suggested and some of these reactions 
with ^^Ca beam on targets of ^^^Th, ^^^C/, 242,244 '^^^Cm nuclei are recently 
used successfully at the JINR Dubna cyclotron U-400. According to the QMFT, 
all these ^^Ca- induced reactions, resulting in different excited compound systems, 
are cold fusion reactions with an " intermediate" amount of excitation energy, com- 
pared to other cold fusion reactions based on Ph beams or light nuclei beams such 
as Mg, Si, and S. The main result of the QMFT is that these three types of 
reactions with different excitation energies refer to three different minima in a 
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potential energy surface of the same excited compound system having different 
depths, and are due to the shell effects of one or both the reaction partners. The 
re-calculated potential energy surfaces using recent binding energy data re-ensure 
the 1977 result. In fact, the QMFT favours the use of ^^Ca beam for these heavy 
isotopes and ^^Ca for the lighter isotopes like ^^"llO, ^^Hl2, ^^^114 and ^90116. 
However, ^^Ca is a radioactive nucleus and hence at present difficult to handle it 
as a beam. Though the experimental physics of radioactive nuclear beams is quite 
different, and at present difficult, the QMFT treat them at par with stable nuclei. 

We have also studied the a-decay chains of the superheavy nuclei ^''"''110, 
287,288,289-]^ "1^4 and ^^^116, obtained in ground-states after neutron evaporations from 
the above mentioned excited compound systems. The model used is the preformed 
cluster-decay model (PCM) of Gupta, also based on the QMFT. The PCM cal- 
culations are then compared with Dubna experimental data and another recent 
calculation. Both the calculations predict similar trends for ct-decay half-lives 
which is due to their used Q-values. The experimental E'^-values (taken as Qa- 
values) are in general higher and show different trends for atleast some cases and 
are the cause for the missing comparisons between calculations and data by about 
two orders of magnitude. An assuring result of PCM is that a-particle preforma- 
tion factors Po are small, ~ 10~^ to 10~^°, required for such rare decays. The other 
calculation uses Pq = 1- 

Finally, the branching of a-decays to other (theoretically possible) cluster de- 
cays are also studied for the ^^^110 nucleus. Interesting enough, some clusters hke 
49-51 (J ^ have their predicted decay half-lives within the present experimental lim- 
its for nuclei with enough available atoms. Note that here, the not yet observed, 
shell effects of clusters are important. Also, the shell stabilizing effects of deformed 
and weakly deformed Z=108 and 106 nuclei are present in these calculations. This 
means that at present the study of cluster decay of superheavy nuclei is mainly of 
theoretical interest to look for new or already present nuclear structure information 
and perhaps for also pointing out some future possibilities with exotic superheavy 
nuclei. 
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Table 1: The a-decay half-lives {logiQTi{s)) and other characteristic quantities for 

2 

ground-state decays of superheavy nuclei with Z=110-116 calculated on PCM, and 
compared with Dubna data (the a-particle energies and decay times r, taken 
as Q^'^' and Tf , respectively) and GLDM calculations. For PCM, the calculations 

2 

are made at Ra — Ci + C2 — Ct, using binding energies from Audi-Wapstra and 
MoUer et al. An increase or decrease in Ra-veilue does not improve the PCM 



results. 
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1.69x10-^* 
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* This chain needs a further careful experimental investigation [14]. 

t These isotopes are also observed in an independent decay chain of 288114^ shown above. 
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Figure Captions 

Fig. 1. The scattering potential for ct-decay of ^^^116 nucleus, calculated as 
the sum of Coulomb and nuclear proximity potential. The tunnelling path 
for PCM is shown, with the first and second turning points Ra and Rb also 
marked. 

Fig. 2. The fragmentation potentials V{A2), calculated at R — Ct, for various 
excited compound systems with Z=110-116. A2 is the mass of the lighter 
nucleus. The heavier nucleus is shown in parentheses for some cases only. 
The ground-state (g.s.) of the compound system is also marked. 

Fig. 3. The calculated formation yields as a function of the mass of reaction 
partners for the compound systems ^^^116* and '^^'^llA* at temperatures cor- 
responding to Ecm =201.1 and 197.2 MeV. 

Fig. 4. The logarithm of a-decay half-lives calculated on PCM, and compared 
with the GLDM and experimental data (decay times r), plotted as a function 
of the parent nucleus mass for various a-decay chains. The data for all the 
observed decay chains are shown and the chains are identified by numbers. 
For PCM, a calculation for Q"^^^' is also shown for 288,289 j^j^^ ^^^^^ 277j^^Q nuclei. 

Fig. 5. The Q-values used in PCM and GLDM, compared with those measured 
{Ea energies) in a-decay chains, plotted as a function of the parent nucleus 
mass. More than one data plotted means more than one observed chain, 
with the same chain identification as in Fig. 4. 

Fig. 6. The PCM calculated Q-valucs, penetrability P, preformation factors Pq, 
and decay half-lives for some cluster decays of the parents of a-decay chain 
for ^^^110, plotted as a function of the parent nucleus mass. 
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Fig. 1 Kumar et al. The formation and decay of Superlieavy. 
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